Electromagnetic induction and damping 
— quantitative experiments using PC interface 
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A bar magnet, attached to an osciUating system, passes 
through a coil periodically, generating a series of emf pulses. 
A novel method is described for the quantitative verification 
of Faraday's law which eliminates all errors associated with 
angular measurements, thereby revealing delicate features of 
the underlying mechanics. When electromagnetic damping is 
activated by short-circuiting the coil, a distinctly linear decay 
of oscillation amplitude is surprisingly observed. A quanti- 
tative analysis reveals an interesting interplay of the electro- 
magnetic and mechanical time scales. 



I. INTRODUCTION 

Laboratory experiments on Faraday's law of electro- 
magnetic induction most often involve a bar magnet mov- 
ing (falling) through a coil, and studying the induced emf 
pulse. Several parameters can be varied, such as the 
velocity of the magnet, the number of turns in the coil, 
and the strength of the bar magnet. The observed pro- 
portionality of the peak induced emf on the number of 
turns in the coil and the magnet velocity provide a quan- 
titative verification of the Faraday's law. 

Commonly it is found convenient to attach the mag- 
net to an oscillating system, so that it passes through the 
coil periodically, generating a series of emf pulses. This 
allows the peak emf to be easily determined by charg- 
ing a capacitor with the rectified coil output. A simple, 
yet remarkably robust setup which utilizes this concept 
involves a rigid semi-circular frame of aluminum, piv- 
oted at the center (O) of the circle (see Fig. 1). The 
whole frame can oscillate freely in its own plane about 
a horizontal axis passing through O. A rectangular bar 
magnet is mounted at the center of the arc and the arc 
passes through a coil C of suitable area of cross section.^] 
The positions of the weights Wi and W2 can be adjusted 
to bring the mean position of the bar magnet vertically 
below the pivot O, and the position of coil is adjusted so 
that its center coincides with this mean position {0 — 0) 
of the magnet. The angular amplitude can be read by 
means of a scale and pointer. The magnet velocity can 
be controlled by choosing different angular amplitudes. 



'^In our experimental setup in the UG laboratory, the coil has 
a diameter of about 10 cm, about the same length, consists 
of several thousand turns of insulated copper wire, and has a 
resistance of about 1000 n. 




FIG. 1. Experimental details. 

allowing the magnetic flux to be changed at different 
rates. 

It is much more instructive to monitor the induced 
emf in the coil through a PC interface, which can be 
readily realized by low-cost, convenient data- acquisition 
modules available in the market. We have used a mod- 
ule based on the serial interface "COBRA-3" and its ac- 
companying software "Measure" , marketed and manufac- 
tured by PHYWE. ||] We specially found useful the var- 
ious features of "Measure" such as "integrate" , "slope" , 
"extrema" , "zoom" etc. In this article we describe mod- 
ified experiments designed to take advantage of the com- 
puter interface. This allows for a quantitative and peda- 
gogical study of (i) angular (position) dependence of the 
magnetic flux through the coil, (ii) verification of Fara- 
day's law of induction and (iii) electromagnetic damping, 
thereby revealing delicate features of the underlying me- 
chanics. 



II. INDUCED EMF PULSE 

The equation for the induced emf £{t) as a function of 
time t can be written as 



(1) 



expressing the combined dependence on the angular gra- 
dient d^/d9 and the angular velocity 1^(6*) = d9/dt. This 
is reflected in the time dependence of £(t), and a typical 
emf pulse is shown in Fig. 2; the pulse-shape is explained 
below for one quarter cycle of oscillation, starting from 
the extreme position of the bar magnet [9 = 9q). As 
the magnet approaches the coil, the induced emf initially 
rises, then turns over and starts falling as the magnet 



1 




—I 1 1 1 1 1 1 4000 4100 4200 4300 4400 4500 

4000 4100 4200 4300 4400 4500 



FIG. 2. A typical induced emf pulse. 

enters the coil and the magnetic flux begins to saturate, 
and finally changes sign as the magnet crosses the center 
of the coil (6 = 0) where the flux undergoes a maximum. 
Thus £ actually vanishes at the point where the angular 
velocity of the magnet is maximum. 

From Fig. 2 the magnitude of £ is seen to be signifi- 
cant only in a very narrow time-interval of about 200 ms, 
which is much smaller than the oscillation time period 
(T w 2 s). This implies that the magnetic fiux through 
the coil falls off very rapidly as the magnet moves away 
from its center, so that d^/dO is significant only in a very 
narrow angular range (typically about 5°) on either side 
of the mean position. As d^/dd = at 6* = 0, it follows 
that d^/dd is strongly peaked quite close to the mean 
position, which accounts for the emf pulse shape seen in 
Fig. 2. This separation of the electromagnetic and me- 
chanical time scales has interesting consquences on the 
electromagnetic damping, as further discussed in section 
IV. 



III. MAGNETIC FLUX THROUGH THE COIL 

In order to quantitatively study the magnetic flux 
through the coil, the "integrate" feature of the software 
is especially convenient. From Eq. (1) the time-integral 
of the induced emf directly yields the change in magnetic 
flux A<i> corresponding to the limits of integration. If 
the lower limit of integration U corresponds to the ex- 
treme position of the magnet {9{ti) = 6*0), @ where the 
magnetic flux through the coil is negligible (valid only 
for large 6*0), the magnetic flux <f>(6') for different angular 
positions 9{t) of the magnet is obtained as 

'Pit) ^ £it')dt' . (2) 

Figure 3 shows a plot of <I>(t) vs. t for a large angular 
amplitude (^o ~ 30°). The time interval during which 
^(t) is signiflcant {^-^ 200 ms) is a very small fraction of 
the oscillation time period (about 2 sec), confirming that 



FIG. 3. Plot of the magnetic flux <1> through the coil with 
time t, showing the rapid change as the magnet crosses the 
center of the coil. 

the magnetic flux changes very rapidly as the magnet 
crosses the center of the coil. As the angular velocity of 
the magnet is nearly constant in the central region, the 
time scale in Fig. 3 can be easily converted to the angular 
{6 = ujt) or linear {x — R9) scale. The points of inflec- 
tion, where d^/dt (and therefore d^/d0) are extremum, 
are at 4200 and 4300 ms, precisely where the peaks occur 
in the emf pulse (Fig. 2). 

IV. VERIFICATION OF FARADAY'S LAW 

For 6*0 >> 5°, the angular velocity of the bar mag- 
net is very nearly constant in the narrow angular range 
near the mean position, and hence the peak emf £niax is 
approximately given by 

\dO J n,ax 

The maximum angular velocity Wmax itself depends on 
through the simple relation (see Appendix) 

t^max = ^ sin(6'o/2) , (4) 

where T is the time period of (small) oscillations. There- 
fore if 9o/2 (in radians) is small compared to 1, then 
<^inax is nearly proportional to 0o, and hence £max ap- 
proximately measures the angular amplitude ^o- 

Eq. (3) provides a simple way for students to quan- 
titatively verify Faraday's law. A plot of the peak emf 
£miix (conveniently obtained using the software feature 
"extrema") vs. Wmax (evaluated from Eq. (4)) for differ- 
ent angular amplitudes should show a linear dependence 
(for large 6*0). While this behaviour could indeed be eas- 
ily verifled by students, the interesting deviation from 
linearity expected for low angular amplitudes {9o ~ 5°), 
for which the 9 dependence of the angular velocity d9/dt 
is not negligible, turned out to be quite elusive. This 
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FIG. 4. Plot of fmax VS. {d£ / dt)l(^Q, showing the deviar 
tion from a straight hne at low angular amplitudes. 

delicate deviation was presumably washed out by the 
large percentage errors in Oq measurements, especially 
for small angles, precisely where this deviation is more 
pronounced. 

An alternative approach, which eliminates the need for 
measuring the oscillation amplitude ^o, is proposed be- 
low. Taking the time derivative of the induced emf, and 
setting = 0, where the angular velocity is maximum, 
we obtain 



9=0 



(5) 



61=0 



which relates the slope off at the mean position to w^ax- 
As this relation holds for all amplitudes it may be used 
to obtain Wmax for different angular amplitudes without 
the need for any angular measurement. The slope at 
the mean position (near zero crossing) is easily measured 
through linear interpolation (see Fig. 2). Thus, a plot 
of £^max VS. {d£ /dt)g=o should show both features of 
interest — the linear behaviour for large angular ampli- 
tudes and the deviation for very low amplitudes. The key 
advantage of this plot lies in completely eliminating the 
errors associated with measurements of oscillation ampli- 
tudes ^0- 

In this experiment the oscillations were started with 
a large initial angular amplitude, and during the grad- 
ual decay of oscillations, the peak voltages £max (on both 
sides of the mean position) and the slope {d£/dt)o^o were 
measured for a large number of pulses, so as to cover 
the full range from large to very small angular ampli- 
tu des. Fig. 4 shows this plot of the averaged Smax vs. 
a/ {d£ /dt)g=o : clearly showing the deviation from linear- 
ity for low angular amplitudes. To provide an idea of 
the angular amplitude scale, the peak voltage of ~ 6 V 
corresponds to an angular amplitude ~ 35°, so that 
the deviations become pronounced when ~ 5°. 



f—i 1 1 1 1 1 1 1 1 r 



+. <>.. 



J I I I I I I I L 



20 40 



100 120 140 160 180 200 220 
t (sec) 



FIG. 5. The normalized peak voltage Vmax(t)/Vmax(0) vs. 
time t for the short-circuit (-I-) and open-circuit (O) cases, 
showing nearly linear fall off. 



V. ELECTROMAGNETIC DAMPING 

To study the nature of the electromagnetic damping 
in the oscillating system, the coil was short-circuited 
through a low resistance (220 il). The oscillations were 
started with a large initial amplitude (still 9o/2 « 1), 
and the voltage V{t) across the resistor was studied as a 
function of time. As the oscillations decayed, the peak 
voltages Knax for sample pulses were recorded at roughly 
equal intervals on the time axis. This voltage V{t) is pro- 
portional to the current through the circuit, and hence to 
the induced emf £{t). As the peak emf fmax is approxi- 
mately proportional to the oscillation amplitude (except 
when the amplitude becomes too small), a plot of Knax 
vs. t actually exhibits the decay of the oscillation ampli- 
tude with time. 

Although an exponential decay of amplitude is more 
commonly encountered in damped systems, the plot of 
the normalized peak voltage Vmax(i)/^ax(0) vs. t shows 
a distinctly linear decay (Fig. 5). To distinguish the elec- 
tromagnetic damping from other sources (friction, air re- 
sistance etc.) the same experiment was repeated in the 
open-circuit configuration, in which case electromagnetic 
damping is absent. In this case the amplitude decay is, 
as expected, much weaker, but significantly it is still ap- 
proximately linear. 

A quantitative analysis of the energy loss provides an 
explanation for this nearly linear decay in both cases. We 
first consider the electromagnetic energy loss. Neglecting 
radiation losses, the main source of energy loss is Joule 
heating in the coil due to the induced current. Integrat- 
ing over one cycle we have 
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(6) 



where TZ is the coil resistance. This may be further sim- 
plified since d$ /dO is significant only in a narrow angular 
range near 9^0 and rapidly vanishes outside. Now, for 
amplitudes not too small, the angular velocity dO/dt is 
nearly constant (« Wmax) in this narrow angular range, 
and therefore taking it outside the integral, we obtain 



'one cycle ' 



7^ 



(7) 



As the angular integral is nearly independent of the initial 
amplitude and therefore of Wmax, the energy loss per 
cycle is proportional to cj,nax, and therefore to ^/E. On 
a long time scale (i >> T), we therefore have 



(8) 



Integrating this, with initial condition £'(0) = Eq, we 
obtain 



^-y/E oc t 



£° -E 



OC t 

OC t . 



(9) 



indicating linear decay of the peak emf, and therefore of 
the amplitude, with time. 

We now consider the energy loss in the open-circuit 
case, where the damping is due to frictional losses. A 
frictional force proportional to velocity, as due to air re- 
sistance at low velocities, will result in an exponential 
decay of the oscillation amplitude. However, a function 
of the type e^"* did not provide a good fit. On the 
other hand, assuming a constant frictional torque r at 
the pivot, which is not unreasonable considering the con- 
tact forces at the pivot, we obtain for the energy loss in 
one cycle. 
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one cycle 



Td9 = T 400 



OC Ve . 



(10) 



This is similar to the earlier result of Eq. (8) for electro- 
magnetic damping, yielding a linear decay of the oscilla- 
tion amplitude with time, which provides a much better 
fit with the observed data, as seen in Fig. 5. The devi- 
ation from linearity at large times is presumably due to 
a small air resistance term. In fact, if a damping term 
—k'E due to air resistance is included in Eq. (8), the dif- 
ferential equation is easily solved, and the solution pro- 
vides an excellent fit to the data. Finally, another term 
-k"E^I'^ should be included in Eq. (8), arising from 



the reduction in the average centripetal force Mluj^BfJI 
with the oscillation amplitude which decreases the 
frictional force at the pivot due to the reduction in the 
normal reaction. 



VI. SUMMARY 

A pedagogically instructive study of electromagnetic 
induction and damping is made possible by attaching a 
PC interface to a conventional setup for studying Fara- 
day's law. By eliminating all errors associated with an- 
gular measurements, the novel method applied for the 
verification of Faraday's law reveals delicate features as- 
sociated with the underlying mechanics. A quantitative 
analysis of the distinctly linear decay of oscillation am- 
plitude due to electromagnetic damping reveals an inter- 
esting interplay of the electromagnetic and mechanical 
time scales. 



VII. APPENDIX 

If the system is released from rest with an angular dis- 
placement ^Oi then from energy conservation 



Mgl{l - cos 6*0) 



(11) 



where M is the mass of the system, / its moment of 
inertia about the pivot O, and I the distance from O to 
the centre of gravity. For small oscillations, the equation 
of motion is 16 = —{Mgl)6, so that the time period is 
given by 



Mgl 



(12) 



Eliminating I /Mgl from these two equations, we obtain 
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